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, Abstract 

o 

\Q , We use lattice sum rules for the static quark potential to determine the beta-function 

0^ ' for symmetric and asymmetric lattices non-perturbatively. We also study the colour field 

distributions in excited gluonic states. 
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In a lattice calculation, the observables are extracted as dimensionless ratios involving the lattice 
spacing a. It is important to obtain results for a range of (the bare lattice coupling coefficient 
given for SU(N) by 2N/g 2 ) in order to explore the approach to the continuum limit a — > 0. The 
region of small lattice corrections is reached when the lattice spacing a obtained from different 
observables has a common /3-dependence. Usually this has been explored by simulating at two 
different values of (3 and matching the results. Another route to this information is from the 
lattice sum rule approach which allows ||, || a derivative of a lattice observable with respect to j3 
to be evaluated at a fixed /3-value by summing plaquette contributions over all space. Expressing 
this in terms of a lattice spacing a determined from that observable, the /3-function b = dfi/dhia 
can then be extracted 0, Q. Moreover, a generalisation to a lattice with different spacings in 
the four directions coming from a generalised Wilson action of the form fiij^ij gives access || 
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to derivatives with respect to these coefficients /3jj by summing plaquettes in the ij plane. It is 
thus possible to evaluate generalised /3-functions 

S if k = i or j and ^ = U if k ^ i or j (1) 



B In afc 9 In afc 

by working on a symmetric lattice (aj = a) at a fixed (3 value. 

For future comparison, the perturbative series for these quantities in terms of the bare lattice 
coupling a = g 2 /Att = \/tt(3 for SU(2) colour fields are || 

b = 2(S + U) = -0.3715(1 + 0.49a + .. .) U - S = 2(3 f = 2/3(1 - 1.13a + . . .) . (2) 

Here we use lattice sum rule techniques to calculate both b and / non-perturbatively. 

We shall apply this sum rule technique to the potential between a static quark and anti- 
quark at separation R. As well as a sum rule study to obtain the generalised /3-functions from 
the ground state potential, we also explore excited gluonic potentials - both to determine the 
generalised /3-functions and to compare with excited string models of the spatial distribution of 
the colour fields responsible for the gluonic excitation. 



2 Sum rules for potentials 

We consider the static quark potential for separation R which is defined as V(R) in lattice 
units. Then the colour fields associated with this pair of static quarks can be measured using 
plaquettes of appropriate orientation. Sum rules have been derived to relate the sum over all 
spatial positions of these colour fields to V(R) and its derivative ||: 

BV 

V + R— = bY / (£L + 2£ T + 2B T + B L ) (3) 
BV 

V + R M = -4Pf^(£ L -B L ) (4) 
BV 

V-R— = -4(3fY,(£T-B T ) . (5) 

Here L refers to longitudinal and T to transverse with respect to the interquark separation 
axis. The notation is that £ and B refer to the difference of the expectation values of the 
plaquette in the appropriate orientation in the presence of the static quarks and in the vacuum: 
£j = (R\D i0 \R) — (0|Dio|0), etc., with □ = ^Tr(l — [/□). Thus they have the interpretation 
of gauge invariant averages of the fluctuation of squared strengths of the colour electric (£) 
and magnetic (B) fields. Since they are a difference between the value in the presence of the 
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(generalised) Wilson loop representing the static quarks and the vacuum value, either sign is 
possible. 

Because the combination £ ±B corresponds to action/energy, we refer to these sum rules as 
'action', 'longitudinal energy' and 'transverse energy' respectively. 

The sum rules are derived || for torelons where there is no self energy associated with the 
static quarks. For the more practical case of the interquark potential between static sources, a 
self energy term must be included in each of the sum rules. This self energy will be independent 
of the inter-quark separation R and hence can be removed exactly by considering differences of 
the sum rules for two different -R-values. 

The remaining obstacle to an evaluation of the sum rules is the presence of the derivative 
dV/dR. This arises in the sum rule derivation from re-expressing the derivative dV{R)/d(5 at 
fixed number of lattice spacings R. In the limit of small lattice spacing a and large number of 
lattice spacings R, this derivative can be accurately expressed in terms of dV(R)/dR at fixed (3. 
For small R, however, there may be substantial systematic error from relating the /5-dependence 
at fixed R to the iZ-dependence at fixed (3. Thus we prefer to retain accuracy by eliminating 
dV/dR from the sum rules. This leads to the relations which can be used to obtain / and b: 

V{R{) - V(R 2 ) = -2/3/ -B l + £ t - B T ) Rl - ^(£ L - B L + £ T - B t )r 2 (6) 

E(£l + 2£ T + 2B T + B L ) Rl - + 2£ T + 2B T + B l )r 2 ' ( j 

As well as using these sum rules for different values of it! to eliminate the self energy terms, 
it is also possible to use the difference between excited states and ground states. The derivation 
of the sum rules can be extended to the case of an excited state and the expressions are the 
same as those given in Eqs. (^), (Q) and (|5|). The self energy is associated with the static 
quarks themselves and is also the same for excited states and ground states and so cancels 
fully for differences of energies between excited and ground states. In practice we use the E u 
representation of D 4 /j to provide the most accessible excited state on a lattice. For more details 
on these excited potentials see ref [0, g|. 

3 Lattice evaluation 

We use a 16 3 x 32 lattice at f3 = 2.4 for this preliminary study with the SU(2) gauge group. For 
many purposes SU(2) colour provides an excellent test bed for QCD studies. In this work the 
scale as set by the string tension corresponds to a lattice spacing a = 0.6 GeV -1 . Generalised 
Wilson loops of size RxT are evaluated with the spatial paths of length R at t = and T being 
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constructed j|] from fuzzed links (using 2 and 13 recursive iterations of the projected sum of 
4 x straight link + 4 U-bends). We study the plaquette signal for all plaquette orientations at 
time slices t = T/2. The correlation of this plaquette sum with the Wilson loop then gives the 
required squared fluctuation of the colour field. 

For the sum rule study, we need the sum of plaquettes over all space in principle. On the 
lattice, we sum the plaquettes over all transverse space (a 16 x 16 plane) for each tl- The 
sum over tl is then approximated by selecting the region within ±2 lattice units of the Wilson 
loop (ie from = — 2 to R + 2). We checked that this approximation did not introduce any 
significant error since the missing ri region contributes noise and no signal to the correlation. 
The data sample used comes from a comprehensive study of 60 blocks of 125 measurements 
separated by 4 update sweeps (3 over-relaxation plus one heatbath). Error estimates use a full 
bootstrap analysis of these blocks. 

The fuzzed link operator to create a static quark and antiquark at separation R with colour 
field in a state of a given symmetry will have an expansion in terms of the eigenstates of the 
transfer matrix 

|-R)=co|Vo> + ci|^i> + ... (8) 

The required ground state contribution is enhanced as T — > oo but at the expense of increasing 
noise. Using a variational basis of 2 different fuzzing levels, we were able to construct very good 
ground state operators. This is essential since the required ground state correlation (Vq|D|Vo) 
is extracted from lattice observables using (for t = T/2) 

(R p t \R T ) = c 2 e- VoT ((Vo\u\V ) + 2h(V 1 \a\V ) + . . .) (9) 

with 

h = 2l e -(vi-v )t (1Q) 

CO 

We can estimate the excited state coefficient h from the study of the Wilson loop itself since 

W(T) = (Rq\R t ) = c 2 e~ VoT (l + /i 2 + . . .) (11) 

so that h 2 ~ e~ v °( N ~ T ' > W(T)/W(N) — 1, where N is the largest time separation which is 
accurately measured and Vq is the best estimate of the ground state energy. To limit the 
contamination of the plaquette correlation by excited states to less than 10%, we require h < 0.05 
which implies h 2 < 0.0025. For T = 2, we obtain h 2 = 0.0005 for R = 2, for R = 4 we find 
h 2 = 0.0025, while for R = 6 we use three fuzzing levels (40, 16 and 0) as basis to construct a 
ground state with h 2 = 0.0036. Thus in each case, to reduce the excited state contamination to 
~ 10 % , a gap of t = T/2 = 1 is sufficient between the operators and the plaquette insertion. 
For electric plaquettes, their unit extension in the time direction implies that a total length of 
T = 3 is needed. Our results for T > 3 have larger statistical errors than those for T = 3 but 
are completely consistent with them within these errors - so confirming the above analysis. Our 
results are collected in Table 1 (using T > 3). 
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R 



! 2 , R 1 T V(R 2 )-V(R 1 ) 



-b 



f 



2, 1 3 0.1892(1) 

3, 2 3 0.1188(2) 



0.502(9) 

0.357(19) 

0.307(32) 

0.315(64) 

0.362(10) 

0.349(23) 

0.351(25) 



0.64(2) 

0.62(3) 

0.64(8) 

0.64(15) 

0.59(3) 

0.64(7) 

0.61(9) 



4 
5 

4, 2 3 0.2129(3) 
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6, 4 3 0.1629(5) 



Table 1: Generalised /3-functions determined from sum rules. 



The lattice observables in the gauge sector with the Wilson action are known to approach 
the continuum limit as a — > with corrections of order a 2 . These corrections, sometimes called 
lattice artefacts, will make the determination of the lattice spacing a different for different 
physical observables on a lattice. The region of a where such differences are small is called the 
scaling region and we expect (3 = 2.4 to be in that region. Nevertheless, d(3/dlna may give 
different values when different observables are used to define a - and that is one of the goals of 
this study. A particularly severe manifestation arises with static potentials at small separation 
R - since the lattice artefacts are known to behave as a 2 / R 2 and are important for R rs o. 
Indeed we see in Table 1 that there are significant lattice artefacts at R = 1 which provide 
an explanation of the discrepancy in the first row of Table 1. We see consistent results when 
R > 1. This provides evidence that, to our precision, non-perturbative determinations of the 
/3-functions are independent of observable for potentials with R > 1 at (3 = 2.4. 

The bare value of a is 0.13, which gives next to leading order perturbative values of b = 
—0.396 and / = 0.85. The effective coupling is expected to be approximately twice as big 
which will decrease the perturbative estimate for / and improve the agreement with our non- 
perturbative result. For b, however, such an increase in a will make the agreement even worse. 
Thus a perturbative evaluation of b is unreliable. 

The /3-function has also been studied non-perturbatively on a lattice by matching the mea- 
sured potentials at two (3 values. Between j3 = 2.4 and 2.5, this gave Q b = —0.277(4). The 
value determined here is the derivative at (3 = 2.4 rather than coming from a finite difference. 
The qualitative features of the two approaches are in agreement, however. 

A non-perturbative study of scaling [|l0| from a thermodynamic approach gave values of 
b = —0.30 and / = 0.66 at f3 = 2.4, although systematic errors are not easy to quantify since an 
ansatz for the functional dependence was assumed. 

We conclude that the generalised /^-functions obtained from sum rules using differences of 
interquark potentials with R > 1 at (3 = 2.4 are consistent with each other, so confirming that 
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these observables have a common a-dependence. The values we obtain are b = —0.35(2) and 
/ = 0.61(3). This a-dependence is not that obtained by the first two terms in the perturbative 
expression. 



4 Excited gluonic potentials 



We measure excited potentials for static sources in the E n representation of the group D 4 ^ which 
is the lattice symmetry |?J of the potential with separation R along a lattice axis. We used 
operators which are U-shaped (ie with links joining the sources of shape n — U) made of fuzzed 
links(with 13 iterations) with transverse extent one and two lattice units. A variational method 
was used to find the optimum combination for the lowest energy state in the E u representation. 
We find the lowest contamination (h? = 0.018 at T = 2) for R = 4. This implies that the study 
of the colour flux distribution in this case will have a contamination of order 25% if T = 3 is 
used in the plaquette correlation approach as discussed above. 

We consider the difference of the symmetric potential (Ai s representation of D4/J at R = 2 
and this E n potential at R = 4. This difference will cause the self-energy terms to cancel. The 
result for this potential difference in lattice units is Ve u (4) — (2) = 0.754(2) yielding, with 
T = 3: b = —0.54(6) and / = 0.68(5). Results for T > 3 are consistent with those for T = 3 but 
have very large errors except for the action sum rule where we do see a significant discrepancy. 
These results for the generalised /^-functions are less precise than those obtained from differences 
among the A\ g potentials (see Table 1) and are not fully consistent. This could arise in part from 
larger lattice artefacts for this particular observable which would be an interesting conclusion. 
This discrepancy can be attributed, however, to the larger contamination of excited states in 
the E u case (as discussed above). One measure of the size of these excited state contributions 
is to compare the T = 2 result with T = 3 and we find much larger differences for E u than for 
Aio. From the T = 4 data for the action sum rule (which controls the value of b) for the E n 
case, we also conclude that it has not achieved a plateau by T = 3. 

Since we have measured the plaquette expectation values in the gluonic excited state (E u ), we 
can explore the spatial distribution of the colour field and make comparison with models for the 
gluonic excitation such as string and flux-tube models. Such comparisons have been made for the 



ground state (Ai ff ) already [jllj, Comparison with such models is most appropriate at large 
R where string-like behaviour is expected. Here we are able to reach R = 8 which corresponds to 
a separation of order 1 fm. Unfortunately, the contamination from excited states (here we mean 
higher energy E n representation states) increases rapidly with R. As an exploratory study, we 
report our results using T = 3 for the generalised Wilson loop since this is the largest T-value 
with a reasonable signal. We specialise to the midpoint (tl = R/2) to minimise the effect of 
self-energy components. The colour field distributions may then be interpreted as applying to a 
state which is predominantly the required lowest energy E n state but with some contamination 
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from excited states. Results for the transverse dependence of the colour flux at R = 8 are shown 
in Figure 1 where they are compared with the symmetric (Ai ff ) potential case. A qualitatively 
similar behaviour is seen for R = 6. 



The interpretation of these results is facilitated by the sum rules. Consider an interquark 
potential given by V(R) = e/R + KR + Vo, where the coefficient e « —0.25 for the A\ g potential 
(this is the 'Coulomb' coefficient) while in string models we expect e = ir for the first excited 
string mode (which corresponds § to the E u case). Then the 'transverse energy' sum rule 
[Eq. (5)] will relate the sum over transverse energy fluctuations to e. This explains the much 
larger transverse energy seen in Figure 1 for the excited gluonic case. The other two sum 
rules suggest that the action and longitudinal energy will have fairly similar contributions when 
summed over rj- for the two cases (the left hand side is 2K in each case) . In fact we see a flatter 
tt distribution for the E u case but the integrals over tt are more alike. This flatter distribution 
in ry corresponds to a 'fatter' flux tube for the excited gluonic case. A more detailed discussion 



will be found in [14] 



In a string model, such as the Isgur-Paton flux tube model |l3|| , the excited modes of the 
string (which contribute to the E u potential) have wavefunctions which are suppressed at = 0. 
This gives a 'fatter' distribution. Unfortunately, in the simplest such picture, this applies only to 
the transverse energy while the other two combinations (action and longitudinal energy) would 
be little changed. 

In conclusion, lattice techniques are able to explore the structure of string-like states in both 
the ground state and excited states. This gives support for such a string- like picture in general 
terms but models are not yet successful at reproducing our results in detail. 

The authors wish to acknowledge that these calculations were carried out at the Centre for 
Scientific Computing's C94 in Helsinki and the RAL(UK) CRAY Y-MP and J90. This work 
is part of the EC Programme "Human Capital and Mobility" - project number ERB-CHRX- 
CT92-0051. 
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Figure 1: The colour flux contributions corresponding to the action (S), longitudinal (El, 
left plot) and transverse energy (Et, right plot) sum rules of Eqs. (3-5) for the static quark 
potential. These are shown in lattice units (with a ~ 0.6 GeV -1 ) versus transverse distance r 
at the mid-point (jl = -R/2) for separation R = 8. The data are for the symmetric ground state 
(Ai 5 representation) and first gluonic excitation (E u representation). 
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